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Abstract. We propose two statistical methods to derive oxygen abundances in HII regions and star forming galaxies and 
calibrate them with a sample of several hundred giant HII regions in spiral and blue compact galaxies as well as of galaxies 
from the Sloan Digital Sky Survey. We show the advantages of our new abundance indicators over previous ones. 
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1. Introduction 

Emission lines due to photoionization of nebulae by massive 
stars are the most powerful indicators of the chemical evolution 
of galaxies in both the near and the intermediate redshift uni- 
verse. The wealth of new data coming from large/deep galaxy 
spectroscopic surveys and from observations with integral field 
units on nearby galaxies should allow a detailed inventory of 
the chemical composition of star forming galaxies. 

Abundance determinations in HII regions are relatively 
straightforward if the electron temperature, T e , can be mea- 
sured directly from the observations (although strong biases 
may affect abundance derivations in high metallicity objects, 
even if lines used for T e measurements are observed, see 
Stasiriska 2005). In many cases, however, T e diagnostic lines 
are too faint to be detected. Pagel et al. (1979) and Alloin et al. 
(1979) pioneered methods that allow one to estimate the metal- 
licity 1 using strong lines only. These methods assume that all 
HII regions are essentially characterized by their metallicity, 
other parameters such as the hardness of the ionizing radiation 
and the ionization parameter being actually linked to metal- 
licity. Empirically, this proves true, at least statistically, as has 
been shown by many theoretical and observational studies. The 
biggest problem is the translation of an observed line ratio into 
metallicity. We will not review here the large amount of litera- 
ture devoted to this subject. References can be found in Perez- 
Montero & Diaz (2005), Pettini & Pagel (2004, hereafter PP04) 
and Stasiriska (2004). There are two important issues. One is 
the choice of the metallicity indicator, the other is the proper 
abundance calibration. Here, we mostly discuss the first of the 
two, as we believe that time is not ripe yet to discuss the sec- 
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ond: one really must understand the physics of metal-rich HII 
regions to be confident in any calibration at high metallicity. 
This requires the systematic use of mid-infrared data to effi- 
ciently constrain the modelling of metal-rich HII regions. 

The most widely used indexes are the following: O23, 
i.e. ([Om] /14959,5007+[Oii] ^3727)/HyS, introduced by Pagel 
et al. (1979), then refined by Mc Gaugh (1991) and 
Pilyugin (2001a); 3 N 2 , i.e. [Ora] ^5007/[Nn] ^6584 intro- 
duced by Alloin et al. 1979 and reconsidered by PP04; 
N2, i.e. [Nn] /16584/Hff first proposed by Storchi-Bergmann 
et al. (1994) and Van Zee et al. (1998), and last recali- 
brated by PP04. Another index that is becoming popular, al- 
though it requires observations in the far red, is S 23 , i.e. 
([S in] ,19069+ [S n] A67 1 6,673 1)/Ha (Vilchez & Esteban 1996, 
Perez-Montero & Diaz 2005). The main problem of the O23 
index is that it is double valued. Therefore another index is 
needed to decide whether one has to chose the low-metallicity 
or the high-metallicity solution. This is not practical if deal- 
ing with data sets where both high and low metallicities are 
expected to be present. The S 23 index is better since it is mono- 
tonic up to about solar metallicity, but a bend at higher metallic- 
ities is expected. The O3N2 and N2 indexes, on the other hand, 
have the virtue of being single valued. However, the reasons 
for this are not purely "physical" but are partly due to "as- 
trophysical" causes: it is an observed fact that, as metallicity 
increases, both the hardness of the ionizing radiation and the 
ionization parameter of HII regions tend to decrease, and that 
the N/O ratio increases on average (at least at high metallici- 
ties). The important role of the "astrophysical factors" , espe- 
cially the variation of the N/O ratio should invite us to some 
caution when relying on them to derive metallicities of objects 
that might have a different nature and history than the objects 
used to calibrate the relation. Such might be the case of high 
redshift galaxies as compared to local dwarf galaxies. 
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Fig. 1. a 12+log O/H as a function of log Ar^O^. Red: HII regions from the spiral galaxies sample; blue: HeBCD sample; cyan: 
SDSS DR3 sample. The thick curve represents our calibration defined by Eq. (1). The two thin curves are vertically displaced 
from the thick one by +0.3 dex. b Same for 12+log O/H as a function of log S3O3. The thick curve represents our calibration 
defined by Eq. (2). 



Another problem that affects both the O3N2 and N2 indexes 
is that they use a "low excitation line", [N 11] /16584 which may 
arise not only in bona fide HII regions, but also in the diffuse 
ionized medium. It is estimated that the diffuse ionized medium 
contributes up to one half of the total Ha emission in galaxies 
(Thilker et al. 2002). Therefore, calibrations of the O3N2 and 
N2 indexes based on giant HII regions might not be relevant for 
the interpretation of integrated spectra of galaxies. 

Ideally, one would require from a metallicity index to: i) 
be single valued with respect to metallicity, ii) have a be- 
haviour dominated by a well understood "physical" reason, 
Hi) be unaffected by the presence of diffuse ionized gas, and 
iv) be independent of chemical evolution. We propose two in- 
dexes: A3O3 , i.e. [Arm]/l7135/[Om]/l5007 and S 3 3 , Le - 
[Sm]/l9069/[Om]/l5007 that fulfill these requirements. In 
Sect. 2 we describe the data base sustaining this work, in Sect. 
3, we present our new metallicity calibrators, and in Sect. 4, we 
discuss their merits with respect to the usual ones. 

2. Data samples 

We have built a large - though not exhaustive - data base 
of published reddening-corrected line intensities. Data for HII 
regions in spiral galaxies have been taken from Garnett & 
Kennicutt (1994), Garnett et al. (1997), van Zee et al. (1998), 
Bresolin et al. (1999, 2004, 2005), and Kennicutt et al. (2003). 
We have also used the sample of blue compact galaxies, ob- 
served by Izotov and coworkers for the quest of the pregalac- 
tic helium abundance, and referred to as the HeBCD sam- 
ple in Izotov et al. (2006). The spectra in this sample have 
very good signal-to-noise ratios. Unfortunately, the wavelength 
range does not reach the [S in] /19069 line. Finally, we have 
used data from the SDSS DR3 sample, i.e. data from the 3rd 
data release of the Sloan Digital Sky Survey (Abazajian et al. 
2005) that have been selected by Izotov et al. (2006) and are 
listed in their Table 1 . 

For both the HeBCD and SDSS DR3 samples, the oxy- 
gen abundances were computed using the temperature de- 
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Fig. 2. log A3O3 versus the electron temperature measured 
from the [Om] /14363/5007 ratio. Blue: HeBCD sample; cyan: 
SDSS DR3 sample. The curve shows the logarithm of the 
[Arm] /l7135/[Om] /15007 emissivity ratio at the temperature 
given by the abscissa (the curve has been shifted for easier com- 
parison with observational points). 

rived from [O in] /14363/5007 and following the prescriptions 
given in Izotov et al. (2006). The formal average uncertainty 
in the computed 12+log O/H values is typically 0.09 dex 
for the SDSS DR3 sample and 0.03 for the HeBCD sample. 
Concerning HII regions in spirals, many of which are located 
in the inner parts of the galaxies, only 22 objects out of over 
300 have the intensity of [0 111] /14363 measured. In some cases, 
other weak auroral lines have been measured. We have however 
decided not to adopt any r e -based method, but rather to use the 
strong line method of Pilyugin (2001a) for high metallicities, 
which uses a combination of O23 and [Ora] /15007/[Oii] A3721 
line ratios (the P-method). This method is based on a thorough 
consideration of many factors, it has been widely discussed 
(Pilyugin 2001b, 2003, Cedres et al. 2004) and it is largely 
used. It is not clear at present whether more recent calibrations 
of the method (Pilyugin & Thuan 2005, Pilyugin et al. 2006) 
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Fig. 3. a 12+log O/H vs. log N2 for our data sample. The thick curve represents Eq. (3). b 12+log O/H vs. log O^Ni. The dashed 
line shows the PP04 calibration. The thick line is our proposed calibration (Eq. 4). The thin curves are vertically displaced from 
the thick one by +0.3 dex. 



are more accurate, because of the problem of abundance bias 
at high metallicity discussed in Stasihska (2005). Since, in any 
event, abundances at high metallicities will have to be recon- 
sidered once systematic measurements of mid-infrared lines 
become available, we simply took the most reasonable of the 
widely used methods, i.e. the calibration by Pilyugin (2001a). 
We discarded the objects for which the method gives 12 + log 
O/H < 8.2, which is the limit of validity defined by Pilyugin for 
the high metallicity range. Pilyugin considers his method to be 
accurate within 0.1 dex. However, the accuracy is certainly not 
as good near the bend in the O23 versus 12+log O/H relation. 



3. Ar 3 6> 3 and S3O3 as metallicity calibrators 

We show in Fig. la the values of 12 + log O/H as a function 
of log A3 03. Red dots correspond to HII regions from our spi- 
ral galaxies sample, blue ones correspond to the HeBCD sam- 
ple, and cyan ones to the SDSS DR3 sample. In total, there are 
666 objects in the plot. Apart from a few outliers in the SDSS 
DR3 sample (for which the oxygen abundance has likely been 
grossly underestimated, and which deserve detailed study), all 
the objects from the three samples merge into one monotonic 
sequence, as expected. The Spearman rank correlation coeffi- 
cient, Rs, is 0.81. The observed trend of Ar^Oi, is for a large 
part due to the decrease of T e as metallicity increases, which 
translates into an increase of the [Arm] /T7135 and [Om] /15007 
emissivity ratio (see Fig. 2). The fact that the hardness of the 
radiation field and ionization parameter decrease with increas- 
ing metallicity adds to the observed relation between 12 + log 
O/H and log Arj03, but does not play a dominant role. Also, as 
expected from stellar nucleosynthesis and confirmed by obser- 
vations of metal-poor HII regions (Izotov et al. 2006), Ar and 
O are produced by the same stars, so, unlike O3N2 or N2 , the 
Ar^Oj index is not affected by chemical evolution effects. Note 
that from Fig. 2 one can infer that the largest observed val- 
ues of Ar^Oi seen in Fig. la must correspond either to a rapid 
change in the hardness of the radiation field and/or ionization 
parameter at the high metallicity end, or to a strong temperature 
gradient inside the HII regions. 



From Fig. la, we can define a new metallicity calibration. 
The observed distribution of points in this figure is well fitted 
by the following expression: 

12 + log O/H = 8.91 + 0.34x + 0.27x 2 + 0.20x 3 (1) 

where x- log Ar^Oj,. The heavy line in Fig. la represents this 
curve. The standard deviation between log O/H (O23) and log 
O/H (A3O3) using this curve is 0.23 dex. Of course, Eq. (1) 
may need revision when the physics of metal-rich HII regions 
becomes better understood. 

Fig. lb shows the values of 12 + log O/H as a function 
of log S3 03. Here, the number of data points is much lower 
(283). Especially, there are no data for the HeBCD sample. 
But the Spearman rank correlation coeffcient is still very high 
(Rs =0.84). The trend shown by the observational points is sim- 
ilar to the one seen if Fig. la, with a smaller slope, due to the 
larger dependence of the ratio of [S m] /19069 and [O m] /15007 
emissivities on electron temperature. The observed distribution 
of points can be fitted by: 

12 + log O/H = 8.70 + 0.28x + 0.03x 2 + O.lx 5 (2) 

where x- log S3O3. This is the heavy line in Fig. lb. The stan- 
dard deviation between log O/H (O23) and log O/H (S3O3) us- 
ing this curve is 0.25 dex. 

Both Eqs. (1) and (2) where found by eye estimate. We 
checked that, at least with our sample, the derivations of metal- 
licity using both methods are consistent: there is no specific 
trend in the differences between both values of log O/H and the 
standard deviation of the differences is 0.09 dex (see Fig. 4a). 

4. Discussion 

We now discuss the relative merits of our Ar303 and S3O3 in- 
dexes with respect to O3N2 and N2. Rather than comparing our 
Figs, la and lb with Figs. 1 and 2 of PP04 (who use a dif- 
ferent data base and different derivations of O/H), we plot in 
Figs. 3a and 3b the values of 12 + log O/H vs log N2 and log 
O3N2, respectively, using our samples. Note that our data points 
are more numerous (740 points instead of 137) and better dis- 
tributed in the entire range of metallicities than the data points 
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Fig.4.a 12+logO/H derived from A r 3 3 (Eq. 1) vs. 12+log O/H derived from S 3 3 (Eq. 2). b 1 2 +log O/H derived from A r 3 3 
(Eq. 1) vs. 12+log O/H derived from N2 using the formula of PP04 (Eq. 3). The black lines show the 1:1 relations. 



in PP04. The heavy line in Fig. 3a represents the calibration of 
the N2 index as given by PP04: 

12 + log O/H = 9.37 + 2.03jc + 1.26x 2 + 0.32x 3 (3) 

where x= log A^. As can be seen, their calibration fits our data 
well. On the other hand, their expression relating 12 + log O/H 
with the O3N2 index, shown with the dashed blue line in Fig. 3b, 
is a very bad fit to our data. This stems from the fact that 
their relation is heavily weighted by only 4 objects in which 
the oxygen abundance has been determined from photoioniza- 
tion modelling and was found to be high. Unfortunately, abun- 
dances derived from models are very uncertain if constraints 
are insufficient. It is also very puzzling that those 4 points are 
quite separated from the rest of the data in the (12 + log O/H 
vs. 3 N2) diagram of PP04, while they are not in their (12 + 
log O/H vs. N2) diagram. Anyway, our Fig. 3b clearly confirms 
their conclusion that O3N2 alone is of no use for low metallic - 
ity objects. It can be used for high metallicity objects, though, 
rather than the calibration proposed by PP04 (their Eq. 3) 2 we 
prefer: 

12 + log O/H = 8.55 - 0.25* (4) 

where x= log 3 Ni. However, comparison of Fig. 3b with Figs. 
la,b shows that there is no advantage in using 3 N2- Given the 
other shortcomings of O3N2 (see Sect. 1), it is preferable to use 
our new indexes, especially Ar 3 3 , which does not require a 
large wavelength range. Note that [Arm] was not considered 
as a "strong line" in previous studies but its intensity is com- 
parable to that of [N 11] at low metallicities. The Ar 3 3 index is 
particularly well-suited to study the metallicities of the galaxies 
in the SDSS and at high redshifts. 

The accuracy of the Ar 3 3 , S 3 3 and N2 metallicity indi- 
cators is comparable at low metallicity, and the former two are 
more accurate at high metallicity, as can be judged from the 
distribution of points between the thin curves in Figs, la, lb 
and 2a. But, as mentioned above, for integrated galaxy spec- 
tra, Ar 3 3 and S 3 3 should be superior even at low metallic- 
ity (at least they do not call for a different calibration to take 

2 Note that PP04 use ([Om] i5007H/3)/([N 11] i6584/Ha) while we 
defined 3 N 2 as [Om] /l5007/[Nn] ,16584. 



into account the diffuse ionized medium). On the other hand, 
both Ar 3 3 and S 3 3 demand a reliable reddening correction, 
which is not required for N2. In the case of global spectra of 
galaxies, reddening correction can be done empirically by us- 
ing the observed relation between the Balmer extinction and the 
galaxy excitation class or the discontinuity at 4000A (Stasiriska 
et al. 2004). This avoids the measurements of Ha and H/3 1 in- 
tensities, often somewhat problematic. 

For the high metallicity end, our calibrations of Ar 3 3 , 
S 3 3 rely on the P-method of Pilyugin (2001a). When metal- 
rich HII regions become better understood, a recalibration of 
all the metallicity indexes in a consistent way may be needed. 
However, our discussion of the merits and drawbacks of Ar 3 3 , 
S 3 3 with respect to other methods will remain valid. 
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